New Limits on the Ultra-high Energy Cosmic Neutrino Flux from the ANITA Experiment 
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We report initial results of the first flight of the Antarctic Impulsive Transient Antenna (ANITA-1) 2006- 
2007 Long Duration Balloon flight, which searched for evidence of a diffuse flux of cosmic neutrinos above 
energies of £v — 3 x lO'^ eV. ANITA-1 flew for 35 days looking for radio impulses due to the Askaryan effect 
in neutrino-induced electromagnetic showers within the Antarctic ice sheets. We report here on our initial 
analysis, which was performed as a blind search of the data. No neutrino candidates are seen, with no detected 
physics background. We set model-independent limits based on this result. Upper limits derived from our 
analysis rule out the highest cosmogenic neutrino models. In a background horizontal-polarization channel, we 
also detect six events consistent with radio impulses from ultra-high energy extensive air showers. 



In all standard models for ultra-high energy cosmic ray 
(UHECR) propagation, their range is ultimately limited by 
the opacity of the cosmic microwave background radiation. 
The UHECR energy above which this becomes significant is 
about 6 X lO'^ eV in the current epoch. This cuts off their 
travel beyond distances of order 50 Mpc as first noted by 
Greisen f 1], and Zatseptin and Kuzmin f2j (GZK). As a result 
of this absorption, the UHECR energy above this GZK cutoff 
is ultimately converted to photons, neutrinos, and lower en- 
ergy hadrons. The resulting neutrinos were first described by 
Berezinsky and Zatsepin (BZ) [3] . In standard UHECR source 
models the BZ neutrino fluxes peak at energies about 2 orders 
of magnitude below the GZK energy. Thus a "guaranteed" 
flux of neutrinos at energies of E^, — lO'^^^" eV exists. Its 
detection is one of the clearest ways to reveal the nature and 
cosmic distribution of the UHECR sources pf], which is one 
of the longest-standing problems in high energy astrophysics. 

The ANITA-1 Long Duration Balloon experiment was de- 
signed specifically to search for this cosmogenic BZ neutrino 
flux. ANITA-1 exploits the Askaryan effect, in which strong 
coherent radio emission arises from electromagnetic showers 
in any dielectric medium |5]. The effect was first observed in 
2000 1 6], and has now been clearly confirmed and character- 
ized for ice as the medium, as part of the pre-flight calibration 
of the ANITA-1 payload |7]. A prior flight of a prototype pay- 
load called ANITA-lite in 2003-2004 led to validation of the 



technique and initial neutrino flux limits that ruled out several 
UHE neutrino models [SJ. 

In a previous paper |@], we describe in detail the ANITA-1 
instrument, payload, and flight system. Reference [9] also in- 
cludes details of the instrument performance during the flight, 
estimates of the overall sensitivity of the instrument to neu- 
trino fluxes, and discussions of possible backgrounds. Be- 
cause of the complexity of the flight system and methodology, 
we refer the reader to ref. |9] for more detail when necessary. 

The ANITA-1 payload (Fig. [T]) launched from Williams 
Field, Antarctica near McMurdo station, on December 15, 
2006, and executed more than three circuits of the continent. 
The payload landed on the Antarctic plateau about 300 miles 
from Amundsen-Scott (South Pole) Station, after 35 days 
aloft. Anomalous stratospheric conditions led to a misalign- 
ment of the polar vortex for the 2006-2007 season, and as a 
result the ANITA-1 trajectory spent an unusually large frac- 
tion of the time over West Antarctica where the ice sheet is 
smaller and shallower. In addition, the payload field-of-view 
to the horizon (at a distance of about 650 km at typical alti- 
tudes of 35-37 km above mean sea level, or 33-35 km above 
the ice surface), often included the two largest occupied sta- 
tions in Antarctica, McMurdo and Amundsen-Scott, and thus 
was subject to higher-than-expected levels of anthropogenic 
electromagnetic interference (EMI). Despite these effects, the 
payload accumulated a net exposure livetime of 17.3 days 



2 



with a mean ice depth in the field of view of 1.2 km, com- 
parable to the attenuation length of the ice at sub-GHz radio 
frequencies |[T3l- ANITA- 1 was thus able to synoptically view 
a volume of ice of ^ 1 .6 M km-'. Our volumetric acceptance to 
a diffuse neutrino flux, accounting for the small solid angle of 
acceptance for any given volume element, is several hundred 
km^ water-equivalent steradians at ~ lO'^ eV Jgt]. 




FIG. 1 : View of the ANITA- 1 payload in launch configuration, with 
photovoltaics at the top and bottom, and antenna clusters between. 
The side of each square antenna mouth is about 0.9 m, and the pay- 
load stands about 8 m tall. 



ANITA- I's antennas are 32 dual-linear-polarization, quad- 
ridged horns, each with a field of view which averages about 
50° angular diameter over their 200- 1 200 MHz working band- 
width. The antennas are arranged in upper and lower rings, 
each with 16 antennas at azimuthal intervals of 22.5°. All 
antennas point at 10° below the horizontal, to maximize sen- 
sitivity to the largest portion of the volume near the horizon 
at 6° below the horizontal. The combined view of all anten- 
nas covers the entire lower hemisphere down to nadir angles 
of about 55°, comprising 99.4% of the area within the hori- 
zon. Radio impulses that exceed the ambient thermal noise 
by about 5o in at least four antennas in coincident upper- and 
lower-ring pairs produce a trigger |9], and the entire antenna 
set of waveforms are then digitized and stored for later anal- 
ysis. Thermal noise fluctuations produce random triggers at a 
rate of about 4-5 Hz, yielding a continuous monitor of instru- 
ment health. These events are incoherent in phase and pro- 
duce a completely negligible background to actual coherent 
radio impulses. 

The event analysis is conceptually simple, but requires de- 
tailed calibration of the instrument to achieve good precision. 
In the results reported here, we accepted only events having 
at least six adjacent antennas with detectable signals. The six 
antenna signals are analyzed using a method of pulse-phase 



interferometry to determine the best arrival direction of the ra- 
dio impulse plane wave, and this direction and its associated 
uncertainty is then mapped onto the Antarctic ice surface by 
reference to onboard payload navigation instruments, with an 
angular precision of 0.2° x 0.8° in elevation and azimuth |9]. 



TABLE I: Event totals vs. analysis cuts and estimated signal efficien- 
cies for unblinded ANITA- 1 data set. 



Cut requirement passed: 


total 


Hpol 


Vpol Efficiency 


(0) Hardware-Triggered 


~ 8.2M 








(1) Upcoming plane wave 


32308 


15997 


16311 


0.93 


(2) Impulsive broadband 


19695 


10095 


9600 


0.98 


(3) Isolated from other events 


9 


8 


1 


0.94 


(4) Isolated from camps 


6 


6 





0.96 


(5) Vpol dominant 











0.99 



To minimize bias, the analysis cuts were optimized on a 
10% randomized sample of the entire data set, and the re- 
maining 90% was blinded from the analysts until the cuts were 
fixed. The cuts proceed as follows: (1) Events that do not re- 
construct to a coherent plane wave in arrival direction are re- 
jected as random thermal noise; events that reconstruct from 
directions above the horizon are also rejected. (2) Events that 
reconstruct but have non-impulsive waveforms from relatively 
naiTow-band sources (< 100 MHz) are rejected. (3) Events 
that cluster with one another in source location to within re- 
construction errors projected onto the ice, or 50 km radius, 
whichever is greater, are rejected. True source candidates 
must be single, isolated events. Note that this cut, and the 
"camp cut" that follows it, are largely but not completely re- 
dundant. (4) Events that coincide in source location with any 
known active or inactive station, camp, aircraft flight path, 
or expedition traverse path, to within reconstruction angular 
errors projected onto the ice, or 50 km radius, whichever is 
greater, are rejected as being associated with anthropogenic 
activity. Even inactive camps or those long-abandoned are 
considered a risk, since left-over equipment might serve as a 
site for charge accumulation and associated electromagnetic 
discharges which could be mistaken for signals. (5) Events 
whose radio waveforms are not predominantly vertically po- 
larized (Vpol) are rejected because, from considerations of the 
Askaryan impulse generation process, and the Fresnel trans- 
mission through the ice surface, they cannot originate from a 
particle shower within the ice sheet. Conversely, strongly hor- 
izontally polarized (Hpol) events are likely to originate from 
above the ice from similar considerations. 

Table J] shows the results of the total event sample after 
unblinding, including the signal efficiency for each cut sep- 
arately. The 10% initial sample is included in the totals. Note 
that the isolation cut (3) is the single most stringent criterion in 
rejecting impulsive events, and this shows that the vast major- 
ity of triggers are not single, isolated events. Signal efficiency 
in each case was tested with a simulated event sample injected 
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FIG. 2: Left: Plot of all reconstructed events, in both horizontal and vertical polarization; major Antarctic stations are indicated on the map. 
Right: events remaining after cuts to remove anthropogenic interference. 6 events remain in the horizontally polarized group, but these are 
non-candidates for neutrino events, as discussed in the text. 



randomly into the data stream, and the final energy-averaged 
efficiency of all cuts is estimated to be 81%. 

In Figure |2] we show the before-and-after maps of recon- 
structed ANITA- 1 events superposed on the Antarctic conti- 
nent. The strong correlation to a small number of stations is 
evident. The 6 surviving Hpol events are by contrast widely 
distributed across the continent, with no known camps or 
bases, either current or former, anywhere in their locale. We 
have investigated the possibility of impulsive signals from 
earth-orbiting satellites seen in reflection off the ice surface 
as a source for these events. This hypothesis is ruled out be- 
cause the waveforms for these events do not show any evi- 
dence of differential group delay from ionospheric dispersion 
which is several ns per MHz in the 200-400 MHz frequency 
range where these events have most of their spectral power. In 
fact the signals are all of durations less than 10 ns. We know 
of no other anthropogenic sources for these events. 

With regard to possible physics sources, our simulations 
of the high-frec|uency tail of impulsive geo-synchrotron ra- 
dio emission 1 131 UM UM from ultra-high energy cosmic ray 
extensive air showers (EAS) suggest that these signals may 
be EAS events seen in reflection off the ice surface |@]. Such 
events are expected to be predominantly Hpol because of the 
strong Fresnel reflectivity in the region near Brewster's an- 
gle, and the overall initial preference for Hpol because of the 
more vertical polar magnetic fields. Our simulations predict a 
handful of such events for the flight, all of which arise from 
UHECR EAS with energies above lO'^ eV; however, the un- 
certainties are large |9]. While these events do not constitute 
a background for our neutrino search because of their incor- 
rect polarization, they are a potentially interesting signal in 
their own right. Further analysis, including a search for simi- 
lar events from above the horizon, is in progress. 
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FIG. 3: ANITA- 1 limits based on no surviving candidates for 18 
days of livetime. Other limits are from AMANDA 1 16], RICE 1 17], 
ANITA-lite LSI, Auger L18J, HiRes |20J, FORTE |19]. The BZ 
(GZK) neutrino model range is determined by a variety of mod- 
els EillilllllllliailllS. 



Based on the approach described in Refs. 24], the re- 
sulting model-independent 90% CL limit on neutrino fluxes 
with Standard Model cross-sections |22] is shown in Fig. [3] 
Here we have included the net livetime and 8 1 % analysis ef- 
ficiency. Exclusion of the volume of ice near all camps and 
events reduces the net effective volume by a few percent. We 
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TABLE II: Expected numbers of events A'v from several UHE neu- 
trino models, and confidence level CL = 100(1 — exp(—Ny )) for ex- 
clusion by ANITA- 1 observations. 



Model & references 



predicted A'v CL,% 



Baseline BZ models 

Protheroe & Johnson 1996 [2ri 0.22 19.7 

Engel, Seckel, Stanev 2001 yT] 0.12 11.3 

Barger, Ruber, & Marfatia 2006 [29] 0.38 31.6 

Strong source evolution BZ models 

Engel, Seckel, Stanev 2001 0.39 32.3 

Kalashev fl/. 2002 [^] 1.03 64.3 

Aramo et al. 2005 [261 1.04 64.6 

Barger, Ruber, & Marfatia 2006 [29J 0.89 58.9 

Yuksel & Kisder 2007 [28J 0.56 42.9 

BZ Models that saturate all bounds: 

Kalashev et al. 2002 [23J 10.1 > 99.99 

Aramo et al. 2005 [26] 8.50 > 99.98 

Waxman-Bahcall fluxes: 

Waxman, Bahcall 1999, evolved sources [12] 0.76 53.2 

Waxman, Bahcall 1999, standard [12] 027 23.7 



estimate that experimental systematics such as variations in 
ice properties and calibration uncertainties in the absolute ra- 
dio signal strength lead to uncertainties of order a factor of two 
on the limit. These model-independent Umits are calibrated 
such that a model spectrum that matched the limit over one 
decade of energy would yield approximately 2.3 events; this 
choice is appropriate to smoothly varying models. The limits 
are an average over all three neutrino flavors, as ANITA- 1 had 
roughly equal sensitivity to Ve, V^, Vx, and the flavors should 
be equally mixed to first order via oscillations for all BZ neu- 
trino models. We plot only an approximate set of bands for the 
BZ neutrino models, which are too numerous to individually 
plot here. 

In Table |ll] we give the total number of events expected 
from selected individual ultra-high energy neutrino models 
which are representative of the range of BZ neutrino expecta- 
tions; we also include the case of a model which saturates the 
canonical Waxman-Bahcall flux bounds for both evolved and 



standard UHECR sources UJ]. ANITA- 1 strongly limits the 
highest BZ neutrino models, which require extremely high- 
energy cutoffs in the parent cosmic -ray sources spectral en- 
ergy distribution. Our limits thus suggest that UHECR source 
spectra extending to 10^^ eV are disfavored. ANITA- 1 sensi- 
tivity approaches a class of models here denoted as "strong- 
source evolution" models, which assume that the UHECR 
source evolution follows the cosmic evolution of more en- 
ergetic sources, for example gamma-ray burst host galax- 
ies ll28ll : these mid-range models are constrained at about the 
60% CL but none are ruled out yet. The ANITA- 1 90% CL 
integral flux limit on a pure E^^ spectrum for the energy range 
IQiS-S eV < £v < 1023-5 eV is E^F < 2 x 10"^ GeV cm^^ g-' 
sr-K 

In summary, we have set the strongest bounds to date on the 



ultra-high energy neutrino flux at energies above 3 x lO'^ eV, 
using the radio Cherenkov method via synoptic observa- 
tions of the Antarctic ice sheets from stratospheric altitudes. 
Our methodology appears to have no observed physics back- 
grounds, but may have detected events due to cosmic -ray ex- 
tensive air showers that are easily separated from the neutrino 
events we seek. 
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